function C = jacobian_f_known_landmark(x0, ne_landmark)

% your code here
eps = 1e-3;

for j=1:numel(x0)
   x = x0;
   x(j)=x0(j)+eps;
   y2=f_known_landmark(x, ne_landmark);
   x(j)=x0(j)-eps;
   y1=f_known_landmark(x, ne_landmark);
   
   C(:,j)=(y2-y1)/(2*eps);
end



